DESY 09-095, LTH 833, SFB/CPP-09-56, HU-EP-09/28 



Meson masses and decay constants from unquenched lattice QCD. 

K. Jansen,! C. McNeile,^ C. Michael,^^ and C. Urbach^ 

^ DESY, Zeuthen, Platanenallee 6, D-15738 Zeuthen, Germany^ 
^ Department of Physics and Astronomy, The Kelvin Building, 

University of Glasgow, Glasgow, G12 8QQ, U.I^ 
^ Theoretical Physics Division, Dept of Mathematical Sciences, 
University of Liverpool, Liverpool L69 3BX, U.I^ 
^ Humboldt- Universitdt zu Berlin, Institut fiir Physik 
Mathematisch-Naturwissenschaftliche Fakultdt I, 
Theorie der Elementarteilchen / Phdnomenologie, 
Newtonstr. 15, 124-89 Berlin German^ 

Abstract 

We report results for the masses of the flavour non-singlet light 0^^, 1 , and 1^ mesons 
from unquenched lattice QCD at two lattice spacings. The twisted mass formalism was used with 
two flavours of sea quarks. For the 0++ and 1^ mesons we look for the effect of decays on the 
mass dependence. For the light vector mesons we study the chiral extrapolations of the mass. We 
report results for the leptonic and transverse decay constants of the p meson. We test the mass 
dependence of the KRSF relations. 

PACS numbers: 11.15.Ha , 12.38.Gc, 14.40.Cs 



* Electronic address: karl.jansen@desy.de 



Electronic address: 'c.mcneile@physics.g la.ac.iik| 
t Electronic address: _ cmi@ liverp o ol . ac . ukj 



'Electronic address: Carsten.Urbach@physik.hu-berlin.de 



1 



I. INTRODUCTION 



Modern unquenched lattice QCD calculations include the dynamics of light sea quarks 
(with pion masses below 300 MeV) and use multiple lattice spacings and volumes This 
has allowed calculations of many basic quantities of long lived hadrons that decay via the 
weak force to be computed to high accuracy. Of particular note is that unquenched lattice 
QCD calculations are now making contact with the results of chiral perturbation theory 
calculations [2|, yj, particularly for light pseudoscalar mesons. 

There has been much less work on studying resonances with the latest generation of 
lattice QCD calculations. Some of the most interesting questions in light quark hadron 
spectroscopy are looking for glueball degrees of freedom in the experimental /o mesons and 
looking for experimental evidence for the exotic 1 mesons. There are new experiments, 
such as Gluex j4| and PANDA ^] that will start around 2015, that aim to study hadronic 
resonances. The new hadronic physics experiments will require results from lattice QCD to 



guide their searches 
been reviewed by 



br new hadrons. The lattice results for light resonances have recently 
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10 



ll|. 



In this paper we test basic lattice QCD techniques to study the bi, oq, and p mesons. 
The observation of the decay of the p meson has been a long goal of the lattice community. 
The issue of dealing with the decay of the p meson has stopped many calculations of weak 
decays such as B ^ pue 12|]. In the case of determining | Vub \ from the semi-leptonic 
decay B pue, the simplest thing is to just ignore this decay and focus on B ^ vrz/e. 
However there are some very important reactions such as 5 — > K*'y and B pr/ that have 
no simple equivalent form factors with a meson that is stable under strong decay. The effect 
of the strong decays on these lattice calculations is an unknown systematic error. It is also 
important to understand the effect of strong decay on the p meson for calculations relevant 

to g-2 Q,Q. 

It has been proposed (see 

mm 

for a review) that the ao(980) contains tetraquark 
or molecular degrees of freedom. It is interesting to see whether quark- ant iquark operators 
actually couple to this state in lattice QCD calculations. Understanding whether the ao(980) 
is a tetraquark is important for classifying the /o and mesons into qq or qqqq multiplets . 

First we define some notation. We call the lightest flavour non-singlet states from the 
lattice calculations with J^*" given by 0^^, 1^~, and 1 as the oq, &i, and p mesons 
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FIG. 1: The effect of decays on the energy levels of a resonance at finite volume. 



respectively at the masses used in the lattice calculation. We include the mass of the state 
when we deal with the experimental state, such as ao(980), p(770). 

The plan of the paper is thus. We first discuss some general issues about the effect of 
hadronic decays on mesons. We then describe the details of the lattice QCD calculation 
and report results for the masses in lattice units. In section [IV] we discuss the interpretation 
of the results for the and hi channels. In section |V] we then discuss the results for the 
masses of the vector mesons. We then discuss the leptonic decay constant of the p meson. 
In the penultimate section we test the KRSF relations. In the final section IVIIII we draw 
our conclusions. 



II. GENERIC BACKGROUND TO THE CALCULATION 

At first analysis, it is not clear that the concept of a hadronic resonance makes sense in 
an Euclidean lattice QCD calculation with a finite box size. Naively, the size of the decay 
width could be a measure of the systematic error on the mass of the resonance on the lattice, 

n 

however there are arguments that suggest this is a pessimistic estimate. Michael |15{] reviews 
some of the phenomenology of unstable hadrons and notes that many unstable mesons fit 
well with mesons that are stable under the strong decays, using SU(3) symmetry for example. 
Also Bijnens et al. [l^ obtained acceptable fits to the masses of the light vector mesons with 
an effective theory (but some parameters coming from a model) that didn't include the effect 
of the vector meson decay. 

In figure [1] we show a "picture" of what we expect happens when a resonance {R) decays 
to two mesons Mi and M2 in the lattice calculation. When the mass of the decay channels 
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and the resonance are close there is mixing between them (an avoided level crossing). The 
hadronic decay in figure [1] requires the creation of a quark- anti-quark pair, so it is only 
present in unquenched lattice QCD calculations. 

For an S-wave decay the threshold for decay is Mr = Mi + M2. For a P-wave decay the 
decay products must carry momentum. For example, in the real world the p decays into 
two pions, via a P-wave decay. The threshold for decay at rest is + (^)^ where L i 



IS 



the side of the box, assuming periodic boundary conditions in space. The CERN group 



13, 



18] found excited masses for the p channel that were consistent with 2yml. + (^)^. For 
heavy quark masses it can be more kinematically favourable to study the decay of the p 



meson with one unit o: 
momentum 
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20 



momentum to decay to a pion at rest and a pion with one unit of 
2l|. It may well be that one of the mesons (Mi or M2) in figure [T] is 
also a resonance, in that case there will be second decay. One example of this is one of the 
decays of the bi meson. 



bi — > urn — > n^TTTin) (1) 

Our lattice calculations can in principle test the effect of the opening of decay thresholds, 
because as we lower the sea quark masses in the calculations, the various decays channels 
become open. In practice it may be hard to see the effect of the open decay as the quark 
mass changes, because other systematic errors may change as well. 

Although it appears that S-wave decays are kinematically easier to observe than P-wave 
decays, the oq and bi mesons are noisier than the p meson. The p meson at rest is stable 
to two pion decay in this calculation, so for this state we try to build in the physics of the 
meson decay by studying the chiral extrapolation formulae in section |Vl We also estimate 
the decay transition amplitude directly on the lattice, to gain an understanding of possible 
consequences of the mixing of the p meson with the two pion state. 

The MILC collaboration claimed to see some evidence for the Oq resonance to decay into 



two light hadrons 



22| . Latter work showed that more analysis was required to understand 



the Oo decay in staggered calculations 
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24 



25|. 



26|. The 



Liischer has developed a technique to compute the scattering phase shifts 
method was applied to 2-d theories and the 0^ theory 28 1. We have not investigated 



newer methods 



29 



future. Morningstar 



30| based on Liischer's technique [26|, but plan to do so in the near 



3l| has recently presented a simple example of the basic method in 



quantum mechanics 



32| 



III. DETAILS OF THE LATTICE CALCULATION 



Our lattice calculation uses the twisted mass QCD formalism [33']. Once a single param- 
eter has been tuned, twisted mass QCD has non-perturbative 0(a) improvement 3J]. We 
call this maximally twisted mass QCD (MTMQCD). This 0{a) iraprovement was checked 
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38, 



391], and 



numerically by scaling studies using quenched QCD calculations [35|, 
has recently been checked in lattice perturbation theory [4^. As a prerequisite for large 
scale unquenched calculations, the phase structure of twisted mass QCD has been stud- 



ied 



41 



42 
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4J, |45|]. The twisted mass formalism has recently been reviewed by 
Shindler [46]. 

The ETM collaboration has already published a comparison of the lattice results for m„ 



and against chiral perturbation theory [471, |48|]. Results for the nucleon and A masses 
and a comparison with chiral perturbation theory are reported in {4^. The masses of the 



flavour singlet pseudoscalar mesons have been presented [50| . Light quark masses and d ecay 



5l| 



5311. the form 



constants from a partially quenched analysis have been published from tliis_ data set 
There are ongoing projects to look at the moments of parton distributions 
factor of the pion 5J], and the properties of heavy-light mesons 55|]. For an overview of 
the broad range of physics projects undertaken by the ETM collaboration see the review by 
Urbach 



56|. 



For the gauge fields we use the tree-level Symanzik improved gauge action [5 
includes the plaquette term t/^^^^, and rectangular (1 x 2) Wilson loops U^^^^^ 



which 



(2) 



l<fjL<U 



with bi = —1/12 and bo = 1 — 861. This choice of gauge action was made after a study of 
the phase structure of unquenched QCD with nf=2 mesons. 

The fermionic action for two degenerate flavours of quarks in twisted mass QCD is given 

by 



5-^ = 0^^ X{x) {Dw[U] + mo + «/i75^^)x(a;) 



(3) 
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FIG. 2: Effective mass plot for the charged p correlators (vector coupling) for the Bq ensemble. 
F and L are the fuzzed and local operators respectively. 



with the Pauli matrix acting in the isospin space, /i the bare twisted mass and the massless 
Wilson-Dirac operator given by 



1 nnr 



(4) 



where 



Uj^{x)ip{x+ajl)—'ijj{x) 



and V*'?/'(x) = — 



U^{x—afi)tp{x—afi)—4'{x) 



(5) 

Maximally twisted Wilson quarks are obtained by setting the untwisted quark mass mg to its 
critical value rricr, while the twisted quark mass parameter /i is kept non- vanishing in order 
to work away from the chiral limit. In eq. ([3]) the quark fields x are in the so-called "twisted 
basis" . The "physical basis" is obtained for maximal twist by the simple transformation 



/ ITT \ / ITT 

ij{x) = exp ( -jl^r^ 1 x{x), ip{x) = x{x) exp ( ^75^^ 
In terms of the physical fields the action is given by 



(6) 



St = a'J2 ^(^) (^7,[V, + V;] - ^^,r' (-y V,V; + m„) + /i^ . (7) 



The generation of the gauge configurations is reported in 



48 



58 



59j. The methods 



used to extract the masses and decay constants of the light mesons, from ny=2 unquenched 



twisted mass QCD are described in [47|, |48|. Correlators separated by 10 trajectories were 



used. The ensembles used in this calculation are summarized in table [H We fit a matrix 
of correlators to a factorising fit form [48]. The basis of smearing functions includes local 
and fuzzed operators. The correlators were calculated with all-to-all quark propagators 



computed using the "one-end-trick" |48l. |60|. 



At finite lattice spacing, there is an order a mixing of mesons with different parity in 
MTMQCD. When studying charged mesons this has the consequence that the p and ai 
mesons mix. Assuming we are at maximal twist, the mixing will be of order a, then at 
large t the lightest state, the p meson, will dominate. The p can be created by a vector or 
tensor current so we used a 4 by 4 matrix of correlators (vector /tensor and local/fuzzed). 
We obtain a good fit with one meson state for t/a>l. We checked these fits using either a 
subset of operators or with more states. 

The charged Oq and hi mesons mix under twisting with spin-exotic mesons so we do not 
expect at large t any significant contributions from parity mixing since those states will be 
heavy. For these cases, we fit a 2 by 2 matrix of correlators (local/fuzzed) from t/a > 3 with 
two meson states. 

In table [TTl we report the masses for the Oq, &i and p mesons in lattice units. In figure [2] 
we plot the effective mass plot for the p correlators for the Bq ensemble. 

In section |IV] we process the raw data and convert the results into physical units. To 
convert the results into lattice units we use the scale from the pion decay constant, at 0/3=3,9 
= 0.0855(5) fm and 0/3=4,05 = 0.0667(5) fm. These scales were consistent with those obtained 
from the mass of the nucleon 49 1. 

In table [TTl we also include the lattice masses for the neutral p° operator. In the twisted 
mass formalism the p° and p"*" mesons are not degenerate because of the flavour violation 
from the twisted mass term. The results in table [TTl show that the p° and p~^ are essentially 
degenerate. A theoretical discussion with numerical examples for why this is so, is contained 
in [611. 



As reported in [61|] the main effect of the flavour violation from the twisted mass term 
is in the mass splitting between the mass of the 7r° and 7r+ mesons. This has implications 
for decay thresholds of the p+ and p° mesons. Experimentally the dominant decays of the 
p+ and p^ meson are to Tr+Tr*^ and tt+tt^ respectively. The physical decay of p" to vr^vr" is 
not allowed, because of isospin symmetry, however at non-zero lattice spacing this decay 
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TABLE I: Summary of ensembles used in this calculation. The format of the measurement column 
is number of blocks times block length. 



Ensemble 


/? 




L3 


X T 


Measurements 


1 


3.9 


0.004 


243 


X 


48 


111 X 8 


B2 


3.9 


0.0064 


243 


X 


48 


78 X 32 





3.9 


0.0085 


243 


X 


48 


66 X 32 


Ba 


3.9 


0.01 


243 


X 


48 


38 X 32 




3.9 


0.015 


243 


X 


48 


44 X 32 




3.9 


0.004 


32^ 


X 


64 


81 X 6 


Ci 


4.05 


0.003 


32^ 


X 


64 


64 X 8 


C2 


4.05 


0.006 


32^ 


X 


64 


66 X 8 


C3 


4.05 


0.008 


32^ 


X 


64 


61 X 8 


C4 


4.05 


0.012 


32^ 


X 


64 


40 X 8 



TABLE II: Masses in lattice units for the oq, bi, and p mesons 



Ensemble 




aniao 


amp+ 


am^o 


Bi 


0.702(52) 


0.539(115) 


0.404(22) 


0.391(17) 


B2 


0.685(28) 


0.573(59) 


0.422(9) 


0.434(17) 


Bi 


0.729(24) 


0.619(31) 


0.428(8) 


0.424(14) 


Ba 


0.681(29) 


0.666(34) 


0.438(6) 




B5 


0.746(30) 


0.699(28) 


0.481(7) 




Be 


0.674(29) 


0.636(53) 


0.416(14) 


0.409(21) 


Ci 


0.552(38) 


0.509(45) 


0.335(12) 


0.352(23) 


C2 


0.555(29) 


0.410(29) 


0.337(12) 


0.344(13) 


C3 


0.526(35) 


0.511(26) 


0.345(8) 




Ca 


0.638(32) 


0.545(19) 


0.368(6) 





is allowed in twisted mass lattice QCD. At =3.9 the mass splitting between the vr" and 
7r+ is approximately 50 MeV at /x=0.004 [48]. This should be compared with one unit of 
quantised momentum of 600 MeV and 450 MeV on the 24^ and 32^ lattices respectively at 
p = 3.9. 
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FIG. 3: The ratio of correlators for the (i) connected neutral (ii) connected plus disconnected 
neutral to the connected charged p correlator for the Bi ensemble. T is the tensor and V is the 
vector current. As discussed in section IVII the different currents renormalise differently, which 
explains whether the ratio tends to one or to the ratio of the square of the renormalisation factors. 
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FIG. 4: The ratio of correlators for the (i) connected neutral (ii) connected plus disconnected 
neutral to the connected charged p correlator for the Ci ensemble. The notation is the same as for 
the caption of figure El 



For a study of flavour singlet vector mesons such as the and cj, evaluation of disconnected 
diagrams is required. Earlier lattice work 62| showed that these contributions are small. 
For the tensor coupling of the vector meson, the considerable variance reduction possible 
using MTMQCD has allowed these contributions to be evaluated with some precision for 



9 



the first time 



631] so yielding first principles results on the uj-p mass difference and mixing. 



Here we are discussing the flavour non-singlet mesons. For the neutral p meson there 
are also disconnected diagrams that contribute to the correlators, because the twisted mass 
formalism breaks isospin symmetry at non-zero lattice spacing. These contributions would 
be expected to be small but, to check this, for the Bi and Ci ensembles we computed 
the relevant disconnected diagram for the vector mesons. Because of favourable variance 



reduction 



481 ]. we are able to determine the disconnected contribution rather precisely for 



neutral p correlations using a vector coupling. The results are in figures |3] and HI As we 
explain in section IVIl the neutral and charged vector currents renormalise differently, thus 
explaining that the ratio of correlators tends to something close to 2, rather than 1. The 
neutral and charged tensor current renormalise the same way, so the ratio of correlators is 
close to 1. For both ensembles the disconnected diagrams make a negligible contribution to 
the correlators, so we do not consider their contribution any further. 



IV. RESULTS FOR THE MASSES OF THE oq AND hi MESONS 



The results for the mass of t 



2007 have been reviewed 



m 



le lightest flavour singlet O^"*" meson from lattice QCD up to 



|. The physics goal is to decide whether a qq interpolating 
operator will couple to the experimental ao(980). The basic summary of the older quenched 
work was that qq interpolating operators did not see the ao(980) meson and coupled to the 
higher non- singlet state 



The unquenched calculation by the RBC collaboration [6J] using nf=2 domain wall 
fermions also found a mass close to the mass of the experimental state ao(1450). In update on 
their analysis, that included 5 times the statistics, the RBC collaboration found 1.11(8) GeV 



for the lightest state in the O"*"^ channel [65|]. McNeile and Michael |6^], in an unquenched 
lattice QCD calculation focused on the mass difference (in the hope that systematics cancel), 
between the 1^~ and the 0^^ mesons. Using this mass splitting it was claimed that the 
lightest state in the O^"*" channel was consistent with the ao(980) state. Lang et al. reported 
masses for the lightest flavour non-singlet 0+"^ consistent with the mass of the ao(980) meson, 
from an unquenched lattice QCD calculation using chirally improved fermions [67| . In an 
unquenched lattice QCD calculation with 2+1 flavours of sea quarks, Lin et al. |68| found 
that the lightest ao state to be consistent with the experimental ao(980). 



10 



> 
o 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 



o (3=3.9, L=24 
□ 0^^ (3=4-05, L=32 
o 7i + ri2(3=3.9,L=24 



o 
<> 



o 
<> 



[] 



Q 
<> 



<> 



ao(1450) 



ao(980) 



0.6 



_L 



_L 







0.1 



0.2 0.3 

2 2 

m GeV 

71 



0.4 



0.5 



FIG. 5: Mass of lightest state in 0''"'' channel with the 7rr/2 decay threshold. 

One complication is that experimentally the Oq decays to ni]. In the two flavour world, 
the lightest rj meson is the flavour singlet pseudoscalar meson at the 800 MeV level. It is the 
mixing between the light and strange loops in a lattice calculation that drives the mixing 
between flavour singlet pseudoscalar states r] and r]'. Hence, the decay thresholds will be 
very different for the Uf = 2 and rif = 2+1 calculations that involve decay to a flavour 
singlet pseudoscalar meson. The ETM collaboration has recently published the masses of 



the flavour singlet pseudoscalar meson (called 772) on these ensembles 69j and these results 
will be used to estimate decay thresholds here. In figure Owe plot the Oq data and the decay 
thresholds. 

To learn how to deal with mesons with open decays on the lattice, we need some simple 
test cases to validate the lattice methods. A bad example to study would be the ai(1260) 



because of its large experimental decay width of 250 to 600 MeV [70|]. The 61 (1235) meson 
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FIG. 6: Mass of bi state snd ttuj threshold as a function of square of pion mass. 



is good choice, because most models treat it as a qq state and its width is not too large 
at 142 MeV [70|]. A direct study of the decay transition bi — > utt has been made on the 



lattice with acceptable agreement 



601 ] with the experimental decay width. To illustrate the 



impact of this (S-wave) decay threshold on the bi meson, we can use the pvr decay threshold 
(because the difference between the p and uj masses is shown to be small 63l|). 

In figure [6] we plot our results from the ETM collaboration for the mass of the 6i meson 
with the estimate of the un threshold, as a function of the square of the pion mass. The 
mass of the lightest state in the bi channel is above the decay threshold. This necessitates 
to include the un operators with the bi operators in a variational analysis, which we plan 
to do in future work. 

In figure [7] we plot the mass difference between the mass of the bi and oq meson as a 
function of the square of the pion masses from a collection of recent unquenched lattice 
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FIG. 7: Mass splitting between the bi and oq mesons. T he p lot includes data from the RBC 



collaboration [65], Hadron Spectrum Collaboration (HSC) 
and the ETMC results from this work. 



, the UKQCD collaboration 



661, 



calculations. The fact that the majority of the results show the mass of the oq meson to 
be lighter than the mass of the bi meson is good evidence for the lightest ag on the lattice 
corresponding to the experimental ao(980) state. The Kentucky group have recently stressed 
that the identification of ao(980) state on the lattice requires an understanding of dynamics 
of the strong decay |7l| . 



V. RESULTS FOR THE MASSES OF THE LIGHT l " MESON 

In this section we will discuss the physical results for the mass of the vector mesons. 
There is much more information on effective field theory for the vector mesons, so there is 
more we can do with the chiral extrapolations in the mass of the light quarks. The data for 
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FIG. 8: The mass of the hght vector meson as a function of the square of the light pseudoscalar 
meson. 



the p meson are useful for applications such as the calculation of the vacuum polarization 

n n 

tensor that is part of the QCD corrections to g — 2 |13l . |14I | and the comparison of the 



electromagnetic form factor of the pion with the vector exchange model 54 1 . 

In figure [S] we plot the mass of the lightest vector meson as a function of the square 
of the pion mass. Our lattice data seem high relative to the experimental mass of the p 
meson. A more detailed comparison with experiment requires a discussion of the chiral 
extrapolations. Also the effect of p decay needs to be accounted for. There has been a long 
history of attempts to deal theoretically with the effect of the p decay on the mass of the p 



meson 



Q 72 
fftl 



73|, 



74 



75 



76|. 



In [71] the vector meson mass as a function of the square of the pion mass, was plotted with 



data from lattice QCD calculations that used improved staggered (MILC collaboration |22| ) . 
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FIG. 9: The correlator c^^t) in equation [8] as a function of time. 



and domain wall fermions (RBC-UKQCD [77[). There was reasonable agreement between 
the data from the different formalisms, although the statistical errors need to be reduced on 
some results (including ours). 

Lattice correlators should have a signal to noise ratio which goes like e"^'"*^ for a 



meson of mass ruM 



78|. We have checked that our data at /3 = 3.9 obeys this relation. So 



there is no fundamental problem with the increase in the statistical errors as the mass of the 
light quarks is reduced. On a subset of the configurations we tried a technique called color 
dilution to improve the signal to noise ratio for the connected p correlators [3]. This did 
not reduce the statistical noise. ETMC have used an extrapolation of the partially quenched 



p masses to reduce the statistical errors [80 1. 



At (3 = 3.9 and p = 0.004, we have also estimated the mixing element between and 
TT^TT^ from a correlator ratio using the method described in |20|]. The three point function 
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ratio was computed using 



^ . ,^ ^ (p(0)|7r(t)7r(t)) 

^ (p(0)|p(t))V2(^(0)vr(0)|7r(t)7r(t))V2 ^ ^ 



When the p-mass and vrvr energy are degenerate, for small enough x [20| (where x = {p \ tttt)), 
this ratio can be fitted to the model in eq. M 

C3 (t) ^ xt + const (9) 

The formalism required, where the p-mass and vrvr energy are not degenerate, is discussed 



m 



811 ] ■ The correlator ratio c^it) is plotted in figure [9] for the decay p 7T{k = ^)7r(fc 



— -^). Since the p-mass is somewhat larger (by 0.19 in lattice units) than the lightest two 
pion energy, we plot in the figure a theoretical curve which modifies eq. |9], taking this into 
account, as used in ref. I20I]. This fit to the three point function ratio gives ax = 0.060(15). 
Since on a lattice, energy is not conserved, we have evaluated the transition amplitude to a 
final state with sufficient momentum that its energy is more than that of the p at rest, so 
strictly a zero decay width. So, to compare with experiment, it is optimum to evaluate the 
coupling constant. This may have some dependence on momentum in general, but it is a 
useful point of reference. The g^^^ coupling defined via 



r = (10) 

biT nip 

is found to be 5'p^7r=5.2(1.3). The corresponding value of (^p^vr from the experimental value 
of the p width is 6.0. So we have consistency between the lattice estimate of the coupling 
between p and tttt and that observed. 

Since we measure the strength of the transition from p to tttt on the lattice (namely x), 
we can estimate the mass shift caused by this mixing. Then with a two-state model with 
energy difference A where 



A = ^2 - ^1 = 2^m2 + (^)2 - nip (11) 
with aA = 0.19 in our case, we obtain, using [20|], a shift (downwards for the p) of 



split 



v/AV4 + a;2 - A/2 (12) 



This mixing produces a 4% downward shift in the mass of the p for ensemble Bi using this 
simplified mixing scheme. This shift is comparable to our statistical error for that state. 
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This suggests that the mass of the vector mesons in figure [8] are largely unaffected by the 
two 7r decay. 

This mixing argument can be used to compare expectations between the Bq ensemble 
with L = 32 and that above with L = 24 above. The differences will be that the energy 
gap will become much smaller (A = 0.08) since the minimum momentum is reduced while 
the mixing contribution (x^) will be reduced proportionally to the spatial volume. The net 
effect is a rather similar estimate which is consistent with our results which show that the 
p mass from the Bq ensemble is half-cr higher than for the Bi ensemble. 

We now discuss the chiral extrapolation of the vector masses to the physical point. For 
the case of an effective field theory for vector mesons, the issues in writing down an effective 



field theory are less clear than for pions. A fully relativistic 
vector fields or a heavy meson effective theo ry ( HMET) 



the different effective theories is discussed in 



16 



16 



agrangian can be used for the 



82( 1 . The connection between 



81. 



The most basic effective field theory for the light vector meson predicts that the mass of 
the vector meson depends on the mass of the pion via la, l82] : 



Mp = M° + ciM^ + C2M3 



(13) 



The pions involved in p decay are not so: 



perturbation theory with power counting 16 



't so p — s> nil can not be studied using chiral 
82 1 . However, Bijnens et al. [l^ successfully 



fitted the masses of the light vector mesons p to (p, including electromagnetic effects, using 
HMET but not including the dynamics of the p — tttt decay. 

The p decay will effect the chiral extrapolation model used to extrapolate the mass of the 



p meson. The Adelaide group have studied different regulators [7^, l75|, TsJ for the effective 
field theory of p decay. This produced additional mass dependence at very light pion masses. 

Models for the effect of ttu and tttt contributions to the mass of the p meson have direct 
implications for the mass of the u meson (which has vrp contributions). Hence lattice results 
for the quark dependence of the mass splitting of the a; to p mesons 63| allow further 
constraints to the study of individual terms. 

Bruns and MeiBner jsj] have published a chiral extrapolation formulae for the mass of 



the p meson. The derivation used a modified MS regulator and a power counting scheme. 

ML 



+ c,M: + C2Mt + csM^^ ln( 



(14) 
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The term with the C3 coefficient is due to the self energy (in the infinite volume limit). Bruns 



and Meifiner 



84l | recommend that the size of the q coefficients obtained from the fits to 



the lattice calculations be checked against constraints from low energy effective constants. 
However they only quote, as reasonable, the constraints that | q |< 3. The Adelaide 
group [7^ claimed to know the sign and magnitude of the C2 coefficient (c2 ~ —1.70 GeV 



but Bruns and MeiBner 
Using one loop chira 

method, Hanhart et al. 

GeV. 

Bruns and MeiBner 



84| claim their bounds are more general. 



perturbation theory and a technique called the inverse amplitude 



85 



86] estimate Ci = 0.90 ±0.11 ±0.13 GeV"^ M° = 0.735 ±0.0017 



84| from an analysis of an old lattice QCD calculation by the CP- 



PACS collaboration [87!], found that the curvature from the non-analytic terms can produce 
either an increase or decrease in the vector mass over a simple linear fit. CP-PACS used 
the string tension (440 MeV) to set the lattice spacing [87|, this corresponds to vq ~ 0.54 
fm, roughly 10 % higher than the preferred tq from the pion decay constant. If there is 
any ambiguity in the lattice spacing, then this can hide the curvature from the non-analytic 
terms. 

Unfortunately the size of errors on the p data and the number of points does not allow 
us to include the C2 and C3 coefficients as free parameters. To get some idea of the effect of 
these terms we use the augmented method 88|, |89| where the physics constraints from 
Bruns and MeiBner j84] can be built into the fit with Bayesian techniques. The augmented 

is used to constrain C2 and C3. 



X. 



aug 



i=2 



(15) 



Schindler and Phillips have recently discussed using an augmented to using information 
from effective theories in chiral extrapolations of lattice data. We use the bootstrap method 



to estimate the errors. In principle given the probability distribution^ 



;he errors on the 



parameters can be obtained by integrating the Monte Carlo integrals |89l . |90| . Chen et al. 
checked |91l] that consistent errors were obtained from a bootstrap analysis and from an 
error analysis based on the augmented being a quadratic function of the fit parameters 
around the minimum. 



We also investigated an approach developed by the Adelaide [7^ group. The Adelaide 



method uses a dipole regulator, rather than the MS scheme, to regulate the effective field 
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TABLE III: The p mass from chiral extrapolation from different fit models at /? = 3.9 



Equation 


Model 


nip GeV 


GeV 


MO GeV 


ci GeV-^ 


C2 (Gey) -2 


C3 {GeVy 


[121 


linear 


0.90(4) 


1.13(8) 


0.89(5) 


0.49(26) 






M 


Bruns and MeiBner 


0.90(5) 


1.07(10) 


0.89(6) 


3.5(4.7) 


-0.09(81) 


-0.82(41) 



theory corrections to the p mass T^] • The extrapolation model for the mass of the p meson 
is 

M. = M; + cMl + (16) 



where E,rw and Etttt are the self energies from the vrvr and txuj states. The fit parameters in 
equation [16] are M°, ci and A^^^. The parameter A^^^ is related to A^ro;. The p self energy 
contribution Stttt contains a cut at nip = 2771^^ for the decay p tttt. For the continuum 
integral we used the principle value of the integral when the decay is open. We found that 
our data was too noisy to get stable fits from this method. We were also unable to resolve 
the quadratic C2 term in equation [131 because the error bars were too large. The original 
study |74| of equation [16] used p masses from lattice QCD with 1% errors at a heavier quark 
masses |74| . 

The summary of the final results is in table lllll We use the pion mass of 135 MeV, because 
we don't include any electromagnetism in the lattice calculation. We also extrapolate our 
results to mass of the notional strange-strange pseudoscalar meson (696 MeV). We call this 
the unitary (p analysis. Note that a better approach to the meson within an nj = 2 
formalism would be to treat the strange quark as a (partially quenched) valence quark with 
a sea of light quarks. 

In figure [lO] we plot the linear fit and the extrapolation model in equation [TU 

The lattice data for the vector mesons seem to prefer a smaller lattice spacing than the 



scales obtained from the pion decay constant |47| and the nucleon mass 49 1, this is probably 
because we are missing some of the effect from the p decay and possibly also from the 
dynamical strange quark. 
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FIG. 10: Fit to the mass of the vector meson using a hnear fit in the square of the pion mass and 
equation 1141 at /? = 3.9. Also included in the plot is the first decay threshold to vrvr for L=24. 

VI. THE DECAY CONSTANTS OF THE p AND MESONS 



We first introduce the leptonic decay constant of the vector mesons, such as the p or 



in the continuum 



92( 1 ■ The decay constant of the vector meson V is defined |93| via 



{0\V^\V) = mpfve^ 



(17) 



where the vector current is defined via 



(18) 



There are other possible (slightly different) definitions of the decay constant of the p defini- 
tions, for example as used by Lewis and Woloshyn [92 1. 
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The decay constants of the p and mesons can be extracted from r decay and e 



annihilation (see 



94 



95| for a discussion). 

^expt _ 208 MeV (19) 

/J^* ~ 216(5) MeV (20) 

^e^xpt _ 233 MeV (21) 

The difference between the experimental values of the /^+^* and /^,?^^* is probably due to the 
problems of extracting the parameters of the p meson from experimental data, rather than 
electromagnetic effects that are important for light pseudoscalar mesons Q]. 
The transverse decay constant {fy^p)) of the V meson is defined by 

(0 I ija^yip I V) = ^fv{^^){p^l^u - Puf^ii) (22) 

where a^y = 'i/2[7^,7y]. It is convenient to introduce the tensor current Tj,^ = ipa^^yip. In 
the lattice calculations we do not include any momentum. 

There is no experimental result for the tensor decay constant fyip) for the p or </) mesons. 



However, light cone sum rules require the transverse decay constant of the p meson 97|,l98l] for 



the extraction of |^ from the -B — > p7 and B K*^ decays. The transverse decay constant 



of the p meson is also used in the analysis of other B decays 12| . There have bee n previou s 



94 



99, 



lOG, 



101 | 



lattice QCD calculations of the transverse decay constants of the p meson 

There needs to be a way to estimate the effect of the strong decay of the p meson to two 
vr on the decay constants, in the same way we tried for the p mass in section |Vl A simple 
test is look at the fv decay constant for the p and (f) mesons as these give us an estimate of 
our accuracy. The majority of older lattice QCD calculations concentrated on the ratio of 
to fv. 

There are various correlators that can be used to extract the fv and fy decay constants. 
For example the correlators in equations [231 EH and [251 Our results are based on factorising 
fits to a basis of 4 by 4 smearing functions that include the local operators as matrix elements 
in the smearing matrix, so the operators in equations [231 EH and [251 are included. 

Y,Y.^V,{x,U)V,{^,^)^) - =M (23) 

X fJ,= l 

E 0)t) - 'f'^f^-'-""" (24) 

X /i=l 
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TABLE IV: Summary of the non-perturbative renormalisation factors used in this calculation. 
The C{p) function is the solution, in equation 1301 of the RG equation for the tensor current 



p 


Za 


^T(/i = i) 


Zy 


c(2 GeV) 


3.9 


0.771(4) 


0.769(4) 


0.6104(02) 


1.01 


4.05 


0.785(6) 


0.787(7) 


0.6451(02) 


1.03 



5^J](r,o(x,t.)T,o(0,0)t) 



3mv(/, 



T\2 -mvi^ 



(25) 



X fJ.= l 



The local vector and T^j, tensor currents need to be renormalised. This involves some 
discussion of the twisted mass formalism. We do all our fits in the twisted bases, however 
the identification of states is done in the physical basis 48|. Assuming that the calculations 
are done at maximal twist (see equation [6]), we have 





1 1 

1 " fj. 1 


j ) cont 


= Zy 


1 v^' 


j) twisted lattice 




v 1 


j ) cont 


= Za 




-A pi 1 twisted lattice 


(^1 


rpa 1 
Vfl 1 


j ) cont 


= Zj' 


1 


1 j ) twisted lattice 



(26) 
(27) 
(28) 

where a takes the values of 1 or 2. Given that we found that the disconnected graphs for 
vector mesons were negligible (in section HVT) . then the connected charged and neutral vector 
mesons give us a separate estimate of the decay constants that use different renormalisation 
constants. This is a useful test of the renormalisation and cut off effects. 



The relevant renormalisation factors Zy^ Zt-, and Z a, h ave been computed 1021 . Il03 | 
using the Rome- Southampton non-perturbative rnethod 1041 ]. The Zy factor has also been 
computed using the conserved vector current [48]. It was found that the conserved vector 
current produced a more accurate estimate of Zy than the Rome-Southampton method, 
so we use the result from the conserved current in this analysis. In this paper we use 
the Za and Zt values calculated through the 'p2-window' method without the use of the 
subtraction of O (a^q^) ter ms. In table HVl we summarise the renormalisation factors used in 



this calculation 



102, 



103] 



The value of the tensor current depends on the scale. The tensor current at /x^ is obtained 
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from that at another scale (/i^) by using the renormahsation group equation. 



(29) 



"^^^^^^l+f^l (7.-/3.70) 



TT 



TT 



+ 1 



1 f Oisifj) 



TT 



with 



[(7i - /9i7o)' + 72 + f^ho - Ml - /?27o] 



- =1i 79 = ^ 

li n ^ Pi n 



(30) 



(31) 



Po 

P2 



11 



38 

, 102 Uf 

16 V 3 ^ 



1 /2857 
64 



5033n/ ^ 325n^ 



2 18 54 

The anomalous dimension for the tensor current has been computed by Gracey 



(32) 



105 



1Q6| 



to three loops in the RI' and the MS schemes. 



7o 
7i 

72 



1 

3 

543 - 26n 



/ 



216 

36^2 + 1440(3^/ + 5240nj + 2784C(3) - 52555 



5184 



) 



(33) 



where the value of the standard constant is G = 1.20206. 



For the coupli ng we used RunDec package 



4-loop evolution 



108 



107l | to compute the coupling from Aqcd using 



1091]. There has not been a calculation of the strong coupling using 

(17) (26) MeVfrom 



informat ion f rom these configurations. We used the value of Aqcd = 26 



110[. The QCDSF value is consistent with that from ALPHA 



lll| . that also used 



QCDSF 
Uf = 2 QCD. 

The results for the leptonic decay constant are reported in table |V] and the results for 
the transverse p decay constant are in table IVI[ The decay constants from the neutral and 
charged vector mesons agree within the errors. We now only consider the decay constants 
of charged vector mesons. In figure [11] we plot the decay constant of the vector meson as 
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TABLE V: Summary of the leptonic decay constant of the vector meson for the different ensembles 
from this calculation. 





Charged 


Neutral 


Ensemble 


afv/ZA 


fv MeV 


fv MeV 


afv/Zv 


Bi 


0.13(1) 


234(18) 


252(13) 


0.179(9) 


B2 


0.148(4) 


264(7) 


283(18) 


0.20(1) 


B3 


0.149(4) 


265(7) 


274(14) 


0.19(1) 


Ba 


0.151(4) 


269(7) 






B5 


0.162(4) 


289(11) 






Be 


0.151(7) 


269(12) 


275(19) 


0.19(1) 


Ci 


0.119(10) 


277(24) 


306(22) 


0.16(1) 


C2 


0.117(6) 


272(14) 


291(14) 


0.152(7) 


C3 


0.117(4) 


272(10) 






C4 


0.121(3) 


281(9) 







a function of the square of the pion mass. There is reasonable scaling between the decay 
constants at /? = 3.9 and [3 = 4.05. The data with larger masses also disagree with the 
value of the decay constant of the meson. The has a small decay width (4.26(4) MeV), 
so we might expect to be able get the properties of this meson correctly. However, the is 
considered to be mostly S7^s, so our neglect of the dynamics of the strange quark could be 
important. 

It has been found that chiral perturbation theory is required to extrapolate the decay 
constants of the light pseudoscalar mesons to their values at the physical quark masses 2|, |3| . 
As discussed in section |V] the application of effective Lagrangian techniques to study the 
p meson is problematic because of the large mass of the p meson relative to the chiral 
scale 8j]. There are expressions for quark mass dependence of the vector meson decay 



constants in 



112l |. The corrections due to loops start at m^logm^ and mg^^. Given the size 



of the statistical errors on the decay constants we didn't try to include any chiral corrections 
in the chiral extrapolations. A simple fit, linear in the square of the pion mass, of the /? = 3.9 
data gives fP''y' = 239(18) MeV and f^''^' = 308(29) MeV. 
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TABLE VI: Summary of the transverse decay constant (/y (/^)) of the vector meson. The scale is 
H=2 GeV. 



Ensemble 


0, fJ^ 1 VjT" 

Y 1 


^Tr(2GeV^ MeV 


f^i2GeV) 
fv 




lOSfS^ 


194(15) 


0.83(4) 


Bo 




195(5) 


0.74(2) 




0.111(2) 


198(6) 


0.75(1) 


Ba 


0.113(3) 


203('6) 


0.75(1) 


Rr 


198('^'l 


91 Si's"! 


n 78n 




0.109(5) 


196(9) 


0.73(2) 


Ci 


0.089(7) 


214(18) 


0.77(5) 


C2 


0.081(4) 


193(9) 


0.71(2) 


C3 


0.089(3) 


214(8) 


0.79(2) 


C4 


0.090(3) 


215(7) 


0.76(1) 



At the moment there are no results for the mass dependence of the transverse leptonic 
decay constants from ef fective fie ld theory, however the formalism for tensor currents has 
started to be developed 



113 



114l |. It will be interesting to see the predictions for the mass 
dependence of the ratio of the transverse to leptonic decay constant from effective field 
theory, because this will test whether a chiral extrapolation of the ratio of the leptonic to 
transverse decay constant results in a cancellation of systematic errors as is hoped. 

There has not been a definitive unquenched calculation of the leptonic decay constant of 
the p meson, although there have been many attempts. Lewis and Woloshyn came within 
1% of the experimental result for the fp in a quenched QCD calculation using the D234 
improved action ^2\- Lewis and Woloshyn summarise older quenched calculations 92|. 
SESAM reported leptonic decay constants for vector mesons that agreed with experiment 



at the 20% level from an unquenched lattice QCD calculation with Wilson fermions 115 1. 



CP-PACS [93| from unquenched calculations with the tadpole improved clover action found 
that they couldn't do a reliable continuum extrapolation of fp. CP-PACS [93] found the non- 
perturbative and perturbative renormalisation factors to be very different. QCD SF obtained 
fp = 256(9) MeV from an unquenched calculation with clover fermions 



and Izubuchi 



lOOH. Hashimoto 



65| obtained fp = 210(15) MeV from a = 2 unquenched calculations that 
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FIG. 11: The leptonic decay constant of the vector meson (as defined in equation I17p is plotted 
as a function of the square of the pion mass. The experimental points for the p and (p are also 
included. 

use domain wall fermions. However this calculation also found that r^^^'^ = 0.549(9) fm from 
the mass of the p meson, so we expect that this is the reason for obtaining a number close 
to the physical point. 

In figure [12] we plot the transverse decay constant of the vector meson as a function of 
the pion mass squared in physical units. The ratio of transverse to leptonic decay constant 
is plotted in figure [IS 

A collection of results for the transverse decay constants are presented in table IVIII 
We also present results from using the ratio of tensor to vector correlators in the bootstrap 



analysis, that we call the "ratio method" . In SOj the ETM collaboration presents results for 
in a partially quenched analysis on the same configurations. We see that our result for 
/J (2 GeV) is approximately 30 MeV higher than most previous results. The RBC-UKQCD 
collaboration also report a result for the transverse decay constant of the K* meson. Only 
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FIG. 12: The transverse decay constant of the vector meson is plotted as a function of the square 
of the pion mass. 



QCDSF [100[ | compute /J on its own, all the others compute -j- and then multiply by 
experiment value for fp. 

From what we call the unitary (j) analysis we obtain fj = 170(14) MeV from the ratio 
method and /J = 222(26) MeV from the direct method, both at the scale of 2 GeV. These 
can be compared wi th th e results i n tab le IVIIi 



Cata and Mateu |121| (see also [122| ) have argued that in the large Nc limit that -j- = 
Their result is consistent with the lattice results in table IVIII for both the quenched 
and unquenched results. There is some ambiguity in the large Nc result, because it doesn't 
depend on the renormalisation scale as it should do. There a re als o predictions for the tensor 



decay constants of the excited vector mesons from large Nc 



12l| . that in principle could 



be 



measured in future lattice QCD calculations that use modern variational techniques 123|. 
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FIG. 13: The ratio of transverse to leptonic decay constant of the p meson as a function of the 
square of the pion mass. The transverse decay constant is at the scale of 2 GeV. 

VII. TESTING THE MASS DEPENDENCE OF THE KRSF RELATIONS 



In this paper we have discussed the mass of the p meson, the lepton decay constant 
fp, and the coupling (yfp.^^ for p decay to rni. Perhaps surprisingly there are postulated 



connections between the three constants, that are called the KRSF relationships 



12J,ll25| 



The orig inal deriv ation of the KRSF relations used the application of the PCAC relation to 



p decay 124 . Il25l |. However, the KRSF r elati ons are also predictions of some effective field 
theories of mesons (see Birse for a review 126 |) , suc h as those with "hidden symmetry" 127| 



and the vector realisation of chiral symmetry 128 1. 



Equation [Ml and equation [35] are known as the KRSFl and KRSF2 relationships |127 | 
respectively. 

/p^ = f.g^. (34) 

ml = flgU (35) 



28 



TABLE VII: Summary of results for transverse decay constants of the p and (p meson. We only 
include the result from the finest lattice of Braun et al. 99 1. 



Group 


Method 


/J (2 GeV) 


fl (2 GeV) 


fp 


Ball et al.fllG. 117. 1181 


sum rule 


155(10) 


208(15) 


0.74(3) 


Becirevic et al. [119] 


quenched lattice 


150(5) 


177(2) 


0.72(2)^^ 


Braun et al. [99] 


quenched lattice 


154(5) 


182(2) 


0.74(1) 


QCDSF [120] 


quenched lattice 


149(9) 






QCDSF [100] 


unquenched lattice 


168(3) 






RBC-UKQCD flOll 


unquenched lattice 


143(6) 


175(2) 


0.69(3) 


This work 


unquenched lattice 


184(15) 






This work (ratio method) 


unquenched lattice 


159(8) 




0.76(4) 



We are using the convention where the physical pion decay constant is f^^ = 130.7 MeV. In 
the effective field theory written down by Georgi 1281 ]. there is an additional 2 on the right 
hand size of equation [35] that makes his model not ag ree with experiment very well. The 
KSRF relations can also be analyzed using AdS/CFT 1291 . Il30 ]. 

In figure [14] we plot gpT,,, from equation [3l] and equation [35] for the /? = 3.9 data, as a 
function of the square of pseudoscalar meson. With in the size of the error bars, the value 
of QpTTn is relatively independent of the pseudoscalar mass. This will be a useful test for 
hadronic effective field theories that include the p meson. 

Some of the work on the KRSF relation in effective field theory is u sed as a q ualitative 
guide to building technicolor models of electroweak symmetry breaking 1271 . Il28 ]. 



VIII. CONCLUSIONS 



The first publication from the ETM collaboration showed impressive agreement between 
the predictions of chiral perturbation theory and the lattice results j47|. In this paper we 
have found that getting agreement between the lattice results and the experimental data for 
the p, bi, oq mesons is much harder. The statistical errors on the masses and couplings are 
too large to look for subtle effects in the chiral extrapolation models. More work, using the 
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FIG. 14: Comparing the ^ptttt coupling from the two KRSF relations with experiment and the 
direct determination described in section IVl 

variational basis for the vector, the bi and the states and more statistics will be needed 
to eventually test the various chiral extrapolations. 

We have started to explore using various tools, such as the Adelaide 0, [ts], 13] method, 
computation of decay widths, and looking for avoided level crossings, to study resonance 
mesons on the lattice. Eventually, the issue o f dea ling with resonances in lattice QCD 
will use Liischer's technique [261] and variants of \13l\\ . These methods are computationally 
intensive, so more pragmatic approaches to studying strong decays on the lattice are still 
important at this time. Liischer's technique for resonances was recently applied to the p 
meson by the CP-PACS collaboration [21i] . 
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